Quarkonium as a tool: cold nuclear matter effects 
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Abstract 

We discuss the quarkonium production as a tool for the study of the Quark Gluon Plasma. In particular, we concentrate 
on the Cold Nuclear Matter effetcs. We show that quarkonium production is also useful for the study of Quantum 
Chromodynamics first principles and the nuclear Parton Distribution Functions. 
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1. Introduction 

The interest devoted to quarkonium production has 
not decreased for the last thirty years. It is motivated 
by the search of the transition from hadronic matter to 
a deconfined state of matter, the so-called Quark-Gluon 
Plasma (QGP). The high density of gluons in the QGP 
is expected to hinder the formation of quarkonium sys- 
tems, by a process analogous to Debye screening of the 
electromagnetic field in a plasma Q. 

Nevertheless, the BNL Relativistic Heavy Ion Col- 
lider data on ii+Au collisions [2 1 have also revealed that 
Cold Nuclear Matter (CNM) effects play an essential 
role at these energies. 

All this confirms that, in fact, the interpretation of the 
results obtained in nucleus-nucleus collisions relies on 
a good understanding and a proper subtraction of the 
CNM effects, known to impact the quarkonium produc- 
tion in proton(deuteron)-nucleus collisions where the 
deconfinement can not be reached. 

We will focus on <i+Au data, where no QGP for- 
mation is possible and only CNM effects are in play. 
Two CNM effects have been identified as important for 
quarkonium production: the nuclear absorption, com- 
monly characterized as a final-state effect on the pro- 
duced quarkonium state and parametrised through an 
effective breakup cross section, and the shadowing, i.e. 
the modification of the parton densities in nuclei relative 
to the nucleon, an initial-state effect. 

Moreover, we have recently realized that the theoret- 
ical uncertainty on the shadowing is not limited to the 



uncertainty on the modification of the parton distribu- 
tions: we do not know the specific production kinemat- 
ics at a partonic level, this prevents us to know the ex- 
act momentum fractions of the gluon in the nucleus and 
thus at which xb to evaluate it. 

In this work, we illustrate the importance of cold mat- 
ter effects on J /if/ production in pPb and PbPb collisions 
at the LHC and we recall that a proper subtraction of 
CNM effects on quarkonium production is mandatory 
before any QGP studies at the LHC. 

2. Our approach 

As we have shown in earlier studies 0|4][5j[6), con ~ 
sidering the adequate J '/if/ partonic production mecha- 
nism - either via a 2 — > 1 or a 2 — > 2 process - affects 
the way to compute the nuclear shadowing and thus its 
expected impact on the J/i/j production. In the follow- 
ing, we will refer to the first scenario as the intrinsic 
scheme, and to the second as the extrinsic scheme. Most 
studies on the J /if/ production in hadronic collisions are 
carried out in the intrinsic scheme. They rely on the as- 
sumption that the cc pair is produced by the fusion of 
two gluons carrying some intrinsic transverse momen- 
tum k T . The partonic process being a 2 — » 1 scattering, 
the sum of the gluon intrinsic kj is transferred to the 
cc pair, thus to the J /if/ since the soft hadronisation pro- 
cess does not alter significantly the kinematics. Thus, 
in such approaches, the transverse momentum Pj of the 
J /if/ entirely comes from the intrinsic k T of the initial 
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gluons. However, the average value of kj is not ex- 
pected to go much beyond ~ 1 GeV. So this process 
is not sufficient to describe the Pj spectrum of quarko- 
nia in hadron collisions Q. 

For Pt>2 - 3 GeV, most of the transverse momen- 
tum of the quarkonia should have an extrinsic origin, i.e. 
the J/if/'s Pj would be balanced by the emission of a re- 
coiling particle in the final state. The J/ift would then be 
produced by gluon fusion in a 2 — > 2 partonic process 
with the emission of a hard final-state gluon. This emis- 
sion, which is anyhow mandatory to conserve C-parity, 
has a definite influence on the kinematics of the J/i// pro- 
duction. Indeed, for a given J /if/ momentum (thus for 
fixed rapidity y and P T ), the processes discussed above, 
i.e. the intrinsic g + g — > cc — > J/i//(+X) and the extrin- 
sic g + g — > J I if/ + g, will proceed on the average from 
initial gluons with different Bjorken-jt. Therefore, they 
will be affected by different shadowing corrections. 

In the intrinsic scheme, the measurement of the J/i// 
momentum in pp collisions completely fixes the longi- 
tudinal momentum fraction of the initial partons: 



*1,2 = 



niT 



exp(+y) = x\ 2 (y,P T ), 



(1) 



with m T = ^M 2 + Pj, M being the J/ip mass. 

However, in the extrinsic scheme, the knowledge of 
the y and Pj spectra is not sufficient to determine x\ and 
%i. Actually, the presence of a final-state gluon intro- 
duces further degrees of freedom in the kinematics, al- 
lowing several (x\ , x%) for a given set (y, Pj). The four- 
momentum conservation explicitely results in a more 
complex expression of X2 as a function of (x\,y, Pj). 



X2 = 



x\mj 



■ M 2 



mje?) 



(2) 



Equivalently, a similar expression can be written for 
x\ as a function of (x2,y, Pt)- Even if the kine- 
matics determines the physical phase space, mod- 
els are anyhow mandatory to compute the proper 
weighting of each kinematically allowed (x 1 ,X2). 
This weight is simply the differential cross section 
at the partonic level times the gluon PDFs, i.e. 
g(x l ,/j F )g(x2,pF)do- gg ^j / ^ +g /dydP T dx 1 dx2. In the 
present implementation of our code, we are able to use 
the partonic differential cross section computed from 
any theoretical approach. In this work, we shall use 
the Colour-Singlet Model (CSM) at LO at LHC energy, 
which was shown to be compatible JS] -see Fig. [T] 
for illustration- with the magnitude of the Pj -integrated 
cross-section as given by the PHENIX pp data IfTUll , 
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Figure 1: rfo"* r | ct /rfv|v=o X Br from gg fusion in pp collisions for 
from 200 GeV up to 14 TeV compared to the PHENIX Q0), CDF |Tj] 
and ALICE data. 



the CDF pp data ifTTI and the recent LHC pp data at 
V^aw = 7 TeV. 

To obtain the yield of J /if/ and T in pA and AA col- 
lisions, a shadowing-correction factor has to be applied 
to the J I if/ yield obtained from the simple superposition 
of the equivalent number of pp collisions. This shad- 
owing factor can be expressed in terms of the ratios R A 
of the nuclear Parton Distribution Functions (nPDF) in 
a nucleon belonging to a nucleus A to the PDF in the 
free nucleon: 



ff(x,Q 2 ) 
Afl n ' de " n (x, Q 2 ) 



i = q,q,8 



(3) 



The numerical parameterisation of Rf(x, Q 2 ) is given 
for all parton flavours. Here, we restrict our study to 
gluons since, at high energy, the quarkonia are essen- 
tially produced through gluon fusion Q. Several shad- 
owing parametrisations are available lfT2l IT3l [T4l IT31 . In 
the following, we shall restrict ourselves to EKS98 0131 . 
which is very close to the mean in the current evaluation 
of the uncertainty |fT31 on the gluon nPDF and exhibit a 
moderate antishadowing. 

The second CNM effect that we take into account 
concerns the nuclear absorption. In the framework 
of the probabilistic Glauber model, this effect is usu- 
ally parametrised by introducing an effective absorption 
cross section <x a b s . It reflects the break-up of correlated 
cc pairs due to inelastic scattering with the remaining 
nucleons from the incident cold nucleus. The value 
of <x abs is unkown at the LHC. At high energy, the heavy 
state in the projectile should undergo a coherent scatter- 
ing off the nucleons of the target nucleus 11611 . in con- 
trast to the incoherent, longitudinally ordered scattering 
that takes place at low energies. As argued in lfT7l[T8l . 
this should lead to a decrease of <x abs with increasing 
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^fsNN■ Compilation and systematic study of many ex- 
perimental data indicate that cr abs appears either con- 
stant [19] or decreasing [20 1 with energy. Hence, we 
can consider our estimates J5] [6] of <x abs at RHIC en- 
ergy as upper bounds for the value of <x abs at the LHC. 
We choose three values of the absorption cross section 
that should span this interval (cr abs = 0, 1.5, 2.8 mb). 

3. Results 

In the following, we present our results for the 
J/xjj nuclear modification factor due to CNM effects 
in the extrinsic sheme in pPb and PbPb collisions 
at -\Jsnn = 5.5 TeV: Rab 



dN J A ! *l{N mll )dN ] J* 



where dN^(dN J p '/) is the observed J If yield in AB = 
pPb,PbPb (pp) collisions and (N co n) is the average 
number of nucleon-nucleon collisions occurring in one 
pPb or PbPb collision. Without nuclear effects, Rab 
should equal unity. 
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Figure 2: J /if/ nuclear modification factor versus y in pPb and PbPb 
collisions at ^Jsnn = 5.5 TeV, using EKS98 1 1 3 ] gluon shadow- 
ing parametrisation and three values of cr a i, s (from top to bottom: 
0, 1.5,2.8 mb) in the extrinsic scheme. For PbPb collisions, the y- 
dependence is shown for various centrality selections. 

In Fig. [2] we show R p ?\, versus y. The curve with 
no absorption allows us to highlight the strong rapidity 
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Figure 3: (Color online) J /if/ nuclear modification factor, RpbPb, in 
PbPb collisions at ^snn = 5.5 TeV versus /V pal t, using EKS98 1131 
gluon shadowing parametrisation and three values of the nuclear ab- 
sorption cross section in the extrinsic scheme. RpbPb is shown for 
two different experimental acceptances in rapidity, |v| < 0.9 and 
Me [2.5,4]. 



dependence of the shadowing. We can also notice that 
the shadowing alone should already be responsible for 
a quite large amount of J If suppression, up to 36 % at 
y = 6. This is expected because one accesses at LHC en- 
ergy the region of very small x in the gluon nPDF (down 
to 1(T 5 ). At backward rapidity, we are in the antishad- 
owing region, with the antishadowing peak at y ^ -5. 
Fig.|2]shows that the y-dependence of Rpbpb for the cen- 
trality bins used by the ATLAS Collaboration iGTl is 
similar for all the bins, with a dip at mid-y. This shape is 
the opposite of the one obtained at RHIC energy |5J|5), 
with a peak at mid-y. One can also notice that 7?p b pb de- 
creases at very large rapidity showing that one has gone 
beyond the antishadowing peak. Here, 7?p b p b is system- 
atically smaller at mid-y than at forward-y. This is also 
illustrated on Fig. [3] with the centrality dependence of 
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^PbPb for two regions in y. This behaviour of the CNM 
effects may partially - or completely - compensate the 
opposite effect expected from cc recombination, with 
a maximum enhancement at y = 0. We should how- 
ever note that the recombination of cc pairs is expected 
to be less important if an experimental Pj cut has to 
be imposed because of a limited acceptance at low Pj. 
Overall, one may observe a 7?pbpb rather independent of 
y resulting of two y-dependent effects. 

Recently, we have also paid attention to the prelimi- 
nary data on T production at RHIC. In Fig. [4] we show 
RdAu vs y according to the extrinsic scheme. The results 
are displayed for 3 values of cr a t, s and for three shad- 
owing parametrisations and these are compared to the 
available RHIC data I22ll23l . 
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Figure 4: T nuclear modification factor in dAu collisions at ysAw = 
200 GeV versus y using the following gluon shadowing parametrisa- 
tions: a) EKS98, b) EPS 08, c) nDSg 



4. Conclusions 

In conclusion, we know from lattice QCD that 
strongly interacting matter undergoes a deconfinement 
transition to a new state, the QGR The objective of 
high energy nuclear collisions is to produce and study 
the QGP under controlled conditions in the laboratory. 
The study of quarkonium production and suppression is 
among the most interesting investigations in this field 
since calculations indicate that the QCD binding po- 
tential is screened in the QGP phase. As demonstrated 
in this work, other phenomena will also impact on the 
quarkonium production rates in nuclear collisions at the 
LHC. For instance, nuclear parton distributions modifi- 
cations and CNM absorption of quarkonium must be ac- 
counted for before any QGP studies at the LHC become 
meaningful. E.G.F. thanks Xunta de Galicia and Min- 
isterio de Ciencia of Spain (FPA2008-03961-E/IN2P3), 
N.M. thanks the F.R.S.-FNRS (Belgium). 
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